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Abstract 



Lowest-order cross sections for qq production and annihilation can be ap- 
proximately corrected for higher-order QCD effects by using a corrective K- 
factor. For energies where quark masses cannot be ignored, the if-factor is 
dominated by the wave function distortion arising from the initial- or final- 
state interaction between the quark and the antiquark. We evaluate this K- 
factor for qq production and annihilation in a quark-gluon plasma by taking 
into account the effects of Debye screening through a color- Yukawa potential. 
We present the corrective i^-factor as a function of dimensionless parameters 
which may find applications in other systems involving attractive or repulsive 
Yukawa interactions. Prominent peaks of the iT-factor occur for an attractive 
q-q color- Yukawa interaction with Debye screening lengths of 0.835 and 3.23 
times the Bohr radius, corresponding to two lowest s-wave qq bound states 
moving into the continuum to become qq resonances as the Debye screening 
length decreases. These resonances, especially the cc and the bb resonances, 
may be utilized to search for the quark-gluon plasma by studying the sys- 



tematics of the temperature dependence of heavy-quark pair production just 
above the threshold. 

PACS number(s): 25.75.+r, 24.85.+p, 12.38.Mh, 13.90.+i 
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I. INTRODUCTION AND SUMMARY 



The possible production of a deconfined quark-gluon plasma (QGP) state during high- 
energy heavy-ion collisions [0-0] is expected to be accompanied by Debye screening of the 
color charges of the constituents in the plasma |8|-|T0|,|2|| . The nature and the characteristics of 
the screening depend on the degree of equilibration and the condition of the plasma. Lattice 
QCD calculations also provide information on the extent of color screening [ID]. While the 
actual temperatures that might be attained cannot be conclusively established at the present 
state of the art, reasonable estimates of the temperature and other plasma properties have 
been made for different scenarios and for different experimental conditions |il~ Hl3[l . On the 
basis of these, corresponding estimates can be made for the color screening that would be 
operative for various scenarios. Recent phenomenology of soft-particle production fl4Ul5 



and the "hot glue" partons scenario [ 12f , suggesting a possible initial gluon-rich environment 
in nucleus-nucleus collisions, raise additional questions regarding color screening in gluon- 
rich matter as well. 

The screening phenomenon modifies the interaction between a quark and an antiquark 
placed in the medium and affects their rates of production and annihilation. Many of 
these reactions, such as the electromagnetic annihilation of q-q pairs and the production of 
strange and charm quark pairs from the plasma, provide valuable signals for the presence 
of the plasma It is clearly important to understand the influence of a screened color 

interaction on various cross sections. 

It is well known that the lowest-order qq annihilation and production cross sections need 
to be corrected to take into account initial-state and final-state interactions. At high energies, 



a A'-factor is used to correct the tree- level rates to achieve a match with experiments [16-20 



At relatively low energies near the q-q threshold, the importance of the final-state interaction 
on the production of heavy quarks has been recognized by many authors pl-p6]1. Similar 
effects of initial- and final-state interactions have also been recognized in many other areas 
of physics p^30|. The Gamow factor (in nuclear reactions [p7| , pion production and pion 
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interferometry and pair annihilations [p9]) and the F-factor (in nuclear beta decay p0[ ) 



have been successfully utilized to correct the lowest-order reaction cross sections or decay 
rates. For example, in beta-decay, the probability for the occurrence of low-energy electrons 
is greatly enhanced because of the attractive final-state Coulomb interaction between the 
electron and the nuclear charge, while the probability for the emission of low-energy positrons 
is suppressed by the repulsive final-state Coulomb interaction. The corrective factors can be 
numerically quite different from unity, especially at energies only slightly greater than the 
rest masses of reacting participants. 

We have in an earlier work presented an approximate representation for the i^-factor to 
correct lowest-order qq cross sections to take into account the initial- and final-state color- 



Coulomb interaction without screening At energies where quark masses are important, 
we have determined the i^-factor by the wave function distortion effects resulting from the 
color-Coulomb interaction between the quark and the antiquark. For the relativistic case at 
very high energies, the corrective i^-factor has been determined from perturbative quantum 
chromodynamics (PQCD) |TB|-|T^]. Our approximate representation of the corrective K- 



factor for the case of no screening interpolates between the low- and the high-energy limits, 
following the interpolation procedure of Schwinger [f^ . It has been tested and found to give 



a good description of the E605 data on dilepton production in the Drell-Yan process |33 



and the experimental ratio of a(e + e~ ^hadrons)/a(e + e~ — > fi + ^~) [|34 . 

In this work we evaluate the fC-factor by including the effects of screening arising from the 
consideration of wave function distortion only. As these corrections are most important at 
low energies, we restrict our investigation to non-relativistic equations of motion at present. 
In future work, we hope to examine the relativistic case where the relative energy is much 
greater than the quark masses. 

Attractive and repulsive Yukawa potentials appear in many other areas of physics. To 
make the corrective i^-factor obtained here useful also for reactants in any attractive or 
repulsive Yukawa potential, we present the corrective factor as a function of two dimension- 
less parameters. For attractive screened interactions, we observe interesting peaks of the 



i^-factor as a function of the screening length, which correspond to single-particle qq bound 
states moving into the continuum as the screening length decreases. These peaks occur when 
the Debye screening length X D is equal to 0.835 and 3.23 times the Bohr radius a B . 

The phase shifts and scattering cross sections for the attractive screened Coulomb poten- 
tial have been evaluated earlier by Morse and his collaborator for the scattering of electrons 
from an atom f35|.|36j and also by Calogero [ 37| . By comparing our if-factor results with 



the phase-shifts of Morse and collaborators, we found that the s-wave resonances at E « 
occur at the same ratios of \ D /a B as the peaks of the .fT-factor. Therefore, the occurrence 
of the peaks of the f^-factor and the resonances at E « are the manifestations of the same 
physical phenomenon arising from emergence of the single-particle qq bound states into the 
continuum to become qq resonances. Because they arise from Debye screening in a decon- 
fined quark-gluon plasma at high temperatures, they can be called "screening resonances" 
to be distinguished from other usual qq resonances such as J /if), if}', and T in the confining 
phase. 

The occurrence of the qq screening resonances, especially the cc and the bb screening 
resonances, may provide a way to search for the quark-gluon plasma. Because the screening 
length is inversely proportional to the temperature, the peaks in the i^-factor will appear 
at specific values of the temperature. Therefore, very large enhancement of the production 
of cc or bb pairs near the thresholds will occur at the corresponding temperatures specified 
by the resonance screening length parameter. A systematic study of the dependence of the 
fT-factor for heavy-quark pair production near the threshold as a function of temperature 
will be useful to map out the effect of screening, which will identify important characteristics 
of the deconfined quark-gluon plasma. 

This paper is organized as follows. In Section II, we discuss the phase-amplitude method 
to calculate the i^-factor for a screened Coulomb-type potential. The corrective i^-factor is 
calculated for attractive and repulsive potentials and the results are exhibited in Section III 
as a function of two dimensionless parameters. In Section IV, we calculate the .fT-factors for 
charm and strange pair production in a quark-gluon plasma. We show how the screening 



modifies the lowest-order cross sections for ss and cc production in Section V. Section VI 
summarizes the present investigation and discusses how the qq resonances in the deconfined 
phase may be used to search for the quark-gluon plasma. 



II. THE A'-FACTOR FOR SCREENED COLOR POTENTIAL 



Basic processes involving q-q initial-state or final-state interactions include dilepton pro- 
duction through qq — > 7* — > 11, charm and strange pair production through qq — > g* — ► 
qq (ss, cc) and gg — > qq (ss, cc). 

In the lowest-order reaction cross sections for these processes, the quark and the antiquark 
are described by plane waves, and their mutual initial- or final-state interactions are not 
included. The quark and the antiquark however interact with each other through a color- 
Coulomb strong interaction, with a large coupling constant. While one can obtain the effect 
of these interactions to the next order in perturbation theory, there are situations in which a 
perturbation expansion involving only the lowest two orders may not be sufficient, as in the 
case of large coupling constants and at energies where the quark masses are not negligible. A 
nonperturbative correction, based on the use of the distorted wave function in the presence of 
the color-Coulomb potential, can provide non-perturbative corrections to the cross sections 

If one places a quark and an antiquark in a color medium such as a quark-gluon plasma or 
gluon-rich matter, the interaction between the quark and the antiquark will be screened. In 
the case of a quark-gluon plasma, the degree of screening depends on the temperature. The 
greater the temperature, the greater the density of quarks, antiquarks, and gluons, and the 
greater will be the degree of screening. The effects of Debye screening in the quark-gluon- 
plasma will modify the configuration space representation of the virtual gluon propagator 
between the quark and the antiquark. The color potential connecting the quark and the 
antiquark will be modified from a color-Coulomb interaction into a color- Yukawa interaction. 
(For an introduction, see page 349 of Ref. 0.) The screened interaction will modify the 
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quark- ant iquark relative wave function, and will therefore change the corrective .fT-factor 
which is given by the square of the wave-function at contact. We therefore need to evaluate 
the wave function for a quark and an antiquark in a screened color potential. 
The q-q Schrodinger equation in the relative coordinate r can be written as 

[V 2 + k 2 - 2fiV(r)]ip(r) = , (1) 

where /i is the reduced mass, and the wave number k at infinite separation is related to 
the nonrelativistic kinetic energy E by k 2 = 2fj,E. The screened color potential V(r) which 
describes the interaction between the quark and the antiquark can be written as 0|],[J 

V{r) = -°^- r , (2) 

where the Debye screening length X D is inversely proportional to quark-gluon plasma tem- 
perature T , and the effective coupling constant a e s is related to the strong coupling a s 
by the color factor Cf. 

a eS = C f a s . (3) 

We have included a negative sign on the righthand side of the potential Thus, when 
q and q are in a color-singlet state, Cf is 4/3 [fTSfl , a e fr is positive, and the interaction is 
attractive. When q and q are in a color-octet state, Cf is —1/6 |18| , a e fj is negative, and 
the interaction is repulsive. 

Following Calogero |j37| , we use the phase-amplitude method to solve the Schrodinger 



equation for the s-state. We write the s-wave radial wave function uo(r) = rtp(r) in terms 
of the phase shift function So(r) and the amplitude function ao(r) as 



Uo(r) = ao(r)|cos5o( r ) [krjo(kr)] — sin5 ( r ) [krr] (kr)] > = oto(r) sin(/cr + 5o(r)). (4) 
From the Schrodinger equation (|l]), the equation for the s-wave phase shift function 5o( r ) is 



d5 (r) 
dr 



-jfe-%VXr)sin 2 [Jfer + <f (r)]. (5) 



This is a first-order differential equation which can be solved numerically starting from the 
origin, with the boundary condition 5q(0) = 0. After the function So(r) is determined, the 
amplitude function ao(r) can be obtained from the relation pT?fl 



(0) expj^A;)- 1 jT dr'[2fiV(r')]sm2[kr' + S (r')]\ . (6) 



We then determine the corrective factor K from the amplitude function at large r using 

«o( r = 0) 



K = 

a (r -> ooj 

In terms of the phase shift function So(r), the corrective i^-factor is given by 



(7) 



K = expj-- dr'[2fiV{r')} sm2[kr' + 5 (r')]| . (8) 



It is convenient to express all radial distances in units of the Bohr radius 

1 



a 



/Lt|a e ff| 



(9) 



The Schrodinger equation in terms of x = r/a B contains only the usual 'Coulomb' parameter 

£ = — (io) 

v 

where v = k/fi and £ can be positive or negative, and the dimensionless screening length 
parameter 

V = ^- (11) 

a B 

The i^-factor is therefore only a function of £ and 77; we shall show results for K (£, if) in the 
next section. 

Note that the nonrelativistic energy E is equal to l/(2fia 2 g ^ 2 ), and the wave number k 
is equal to l/(a fl |£|); thus, E and k decrease as |£| increases. In studying K(£,rj) for the 
quark- gluon plasma, one can consider |£| a measure of the inverse of \[E and rj a measure 
of the inverse of the quark-gluon plasma temperature T. 
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III. FEATURES OF THE CORRECTIVE FACTOR FOR SCREENED 

POTENTIALS 



To illustrate some general features of the corrective factor K(£, 77), it is instructive to con- 
sider a simple potential similar to the screened potential (0), for which analytical corrective 
factor can be easily obtained. The truncated color-Coulomb potential V(r) = —6(a—r)a e g/r 
is qualitatively similar to the screened Yukawa potential (Q), and can be chosen to match 
the Yukawa potential (0) with the same volume integral by setting a = \^2\ D . It can be 
shown easily that the corrective i^-factor for this truncated color-Coulomb potential is 



K 



sin(A;a + 6 (a)) 



F Q (£,ka) 



2tt£ 



exp{-2vr£} 



(12) 



where 60(a) is the phase shift for the truncated color-Coulomb potential V(r) and F (£, ka) 
is the usual regular s-wave Coulomb wave function (Eq. (14.1.3) of Ref. |39|] ). 

For large values of a and ka » 1, we have 60(a) = £ ln(2/ca) + argT(l — z£) and 



K 



sin^n 



/ sin 9 + g cos 9 



2vr£ 



l-exp{-27re}' 



(13) 



where 



9 Q = ka + £ ln(2fca) + arg T(l - if) 



(14) 



and the asymptotic expansion of Fo(£, ka) gives [g£ 



f _ 1 j $e-e 

1 2ka 8k 2 a 2 



(15) 



e < 3 ~ 2£ 
2ka 8k 2 a 2 



(16) 



The first factor on the righthand side of Eq. ( |T3"D is close to unity, with small oscillatory 
corrections. Therefore, for large values of a (which is proportional to 77), K(£,r]) approaches 
and oscillates about the Gamow factor 27r£/(l — exp{— 27r£}), provided that ka is also large, 
corresponding to r\ » |£|. 
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For small values of a (or rj) and ka, we can evaluate the phase shift So(a) in the Born 
approximation and obtain 

5 (a) = £JfeV(l - JfeV/6) . (17) 

From Eq. (0), an expansion of the Coulomb wave function in powers of ka leads to 

[1 + £ha{l - k 2 a 2 /Q) - fcV/6] 2 



K 

[1 - £/ca - (2£ z - lj/e z a76j z 

(18) 



1 - 2^2 sign(£)?7 ' 

where sign(£) is the sign of £. By our convention in Eq. fl2|), sign(£) = 1 for attractive 
g-g interactions and sign(£) = — 1 for repulsive q-q interactions. Therefore, for small values 
of 7] and small ka (which corresponds to |£| >> 77), K(£,r}) is approximately unity and is 
independent of |£|. 

Having discussed the limiting behavior of K(£,rj), we return to the screened color- 
Coulomb potential of Eq. (fj) and display the corresponding corrective factor K(£, rj) as 
a function of £ and rj. The corrective factor K(£, rj) for — 1 < £ < 1 and 77 = 0.1, 0.2, 0.7 
is shown in Fig. 1. It is greater than unity for positive £, corresponding to an attractive 
interaction pulling the interacting particles together. It is less than unity for negative £, 
corresponding to a repulsive interaction pushing the interacting particles apart. For small 
77, it approaches unity, indicating the diminishing influence of the interaction between the 
reacting particles for small screening lengths. The values of K(£, rj) for small rj and large £ 
agree with the estimates from Eq. (p~8|) . 

We can also compare our results to the limit of no screening. It may be recalled that in 
this limit with 77 — > 00, the corrective factor K(^,rj) approaches the Gamow factor which 
is only a function of £. The Gamow factor is shown in Fig. 1 as the dashed curve. For the 
range < £ < 1 shown in Fig. 1, K(£, rj) increases as £ and rj increase. It is quite close to 
the Gamow factor for rj > 0.7. 

In Fig. 2 we show K(^, rj) as a function of £ in the range 1 < £ < 3 for various values 
of the screening length parameter 77. The Gamow factor is the dashed curve. For small rj, 

10 



K(^, rf) is nearly constant in £. As the screening length parameter 77 increases, the corrective 
K-f actor increases to become greater than the Gamow factor, reaches a maximum value, 
then reverses itself and oscillates about the Gamow factor. This oscillation about the Gamow 
factor can be understood from the approximate estimate of Eq. (|13|) . 

The peculiar behavior of the rise of K(£, rf) in the region of r\ ~ 1 merits closer scrutiny. 
Accordingly, in Fig. 3a we study K(£,r]) as a function of 77 in the range 0.4 < 77 < 2. We 
find that K(C,,t]) is a maximum at 77 = 0.835. The larger the values of £, the greater is 
the maximum value of K(C,,i]). The peak values of the .fT-factor are much greater than 
the corresponding Gamow factor, as indicated by the ratio -fT(£,?7)/(GaHiow factor) as a 
function of rj in Fig. 3b. 

The prominent peak of K(^, rf) at 77 = 0.835 is due to the emergence of the lowest bound 
state of the system into the continuum to become a qq resonance as the screening length 
decreases. Movement of bound states into the continuum is familiar in the context of a 
single-particle system in a finite-range potential, as in the scattering of an electron from an 
atom P5] , |3"6| ] . An attractive screened potential with a large screening length is able to hold 
many negative-energy bound states. As an illustration, consider an attractive Coulomb 
potential which holds an infinite number of bound states. As the screening length X D 
decreases from infinity, the Coulomb potential becomes a Yukawa potential with a decreasing 
range and its bound states are pushed into the continuum. Using semiclassical estimates, 
it has been noted that when the screening length A is comparable to the Bohr radius, the 
screened potential cannot hold any bound state 00. Numerical calculations show that 
rj « 0.84 is the screening length parameter for which the lowest bound state in the screened 
potential becomes unbound [[|. This coincides with the location of the peaks of K(^,r]) in 
77. Combining the bound-state result of Ref. with our continuum result for the iT-factor, 
we conclude that the peak of the corrective factor at rj = 0.835 arises because the lowest qq 
bound state emerges into the continuum to become a qq resonance as the screening length 
decreases. The occurrence of similar resonances for electron scattering from an atom has 



There are other well-known physical phenomena in which bound states emerge into the 
continuum to become resonances as the range of the attractive potential decreases. For 
example, in nuclear physics the s-wave strength function obtained from very low energy 
neutron scattering shows peak values as a function of the nuclear mass number at A pa 60 
and A pa 160. (See page 230 of Ref. |4T|.) The s-wave strength function is related to the 
wave function at the nuclear surface, which is correlated with the wave function at the 
nuclear interior. The peaks of the s-wave strength function at A pa 60 and 160 correspond 
to the emergence of the bound single-particle 4s 1/2 and 3si/2 states respectively into the 
continuum as the nuclear mass number decreases. (See page 239 of Ref. ||40|| .) The decrease 
of the nuclear mass number is accompanied by a decrease of the range of the attractive 
interaction. As a consequence, nucleon single-particle states are pushed into the continuum 
as the range of the interaction decreases, similar to the emergence of qq bound state into 
the continuum as the screening length decreases, as found in the present study. 

The behavior of the qq resonance at r\ — 0.835 as a function of energy can be studied 
by examining the asymptotic phase shift [So = So(r — > 00)) for the screened potential as 
a function of E (in units of (2/iia^)™ 1 ). For rj > 0.835, there is a bound state and the 
asymptotic phase shift Sq starts at So = it at E = (according to Levinson's Theorem 
and Sq decreases monotonically as the energy increases (as in the short-dashed curve in Fig. 
4). At rj = 0.835, the asymptotic phase shift starts at So = at E = 0, rises to approximately 
7r/2, and then decreases as the energy increases (the solid curve in Fig. 4). When rj < 0.835, 
So rises from zero to a value close to 7r/2 and decreases at higher energies. As 77 decreases, the 
resonance peak (the location of the maximum Sq) moves to higher energies, the maximum 
phase shift decreases, and the width of the structure increases. For 77 < 0.835, although the 
phase shift nearly passes through n/2 the corrective i^-factor does not show a prominent 
peak structure at the resonance energy (the location of maximum So), because the corrective 
.fT-factor is a rapidly decreasing function of the energy (approximately as 1/yE). The best 
indicator for the presence of this qq resonance at 77 = 0.835 appears to be the the prominent 
peak structure of the corrective i^-factor as a function of 77, much as the manifestation of 
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s-wave resonances from the s-wave strength function plotted as a function of the nuclear 
radius or the nuclear mass number. 

It is interesting to inquire whether there are similar peaks of K(£, 77) at other locations 
of 77. We have found that there is another similar peak of K(£,rj) at 77 — 3.23 (Fig. 5) 
corresponding to the second s-wave bound state emerging into the continuum to become 
a qq resonance. It is clear that other peaks of i^-factor can occur at higher values of 77 
when single-particle s-wave bound states are pushed into the continuum to become s-wave 
qq resonances as \ D decreases. 

The .ff -factor for negative values of £ gives the corrective factor for the lowest-order 
cross section involving two particles subject to a mutually repulsive screened interaction, as 
in the case of a q and a q in color-octet states or two 7r + particles subject to a screened 
Coulomb repulsion in pion interferometry. In Fig. 6 we show K(£, 77) as a function of £ in 
the range —3 < £ < for different screening length parameters A D ={0.1, 0.6, .., 3.1}. The 
corresponding Gamow factor (the no-screening limit) is also displayed. As one observes in 
this range of 77, the X-factor for the screened potential is greater than the Gamow factor, 
and the deviation from the Gamow factor becomes greater, the larger is |£|. The screening 
length 77 need to be much greater than |£| to reach the no-screening limit. Thus, the effect 
of screening is also important for repulsive screened potentials, and the use of the Gamow 
factor for the case with screening may sometimes greatly over-correct the initial- or final- 
state interactions if 77 is small. 

IV. CORRECTION FACTORS FOR TYPICAL PLASMA 

In order to study the case of screening in the quark-gluon plasma, we make a minor 
modification to the general case examined in the last two sections. While we use the non- 
relativistic potential (|2|) to describe the interaction between a quark and an antiquark, 
we would like to include relativistic kinematics for the system, as we consider also kinetic 
energies with associated velocities of the order of 0.5. The use of relativistic kinematics 
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will avoid errors due to the neglect of the higher powers of the momentum, which would be 
present if we restricted ourselves to non-relativistic energy-momentum relations. In place of 
Eq. ([!]), the equation of motion to be solved is [f31"1 , |42 



p 2 + 6 2 -2e^oJ^ = (19) 

where Aq(t) = V(r), and the quantities b and e w can be expressed in terms of the center-of- 
mass energy i/s as 

e w = (s - m\ - m 2 2 )/2^~s (20) 

b 2 = [s 2 - 2s(m\ + m 2 2 ) + (mj - m 2 2 ) 2 }/As = el - m 2 w (21) 
where m w = m 1 m 2 /y/s. The corresponding dimensionless parameters £ and rj are 

v b 

and 

V = — = ^oacs^ ■ (23) 
a B 

Values of 1.5 GeV and 0.15 GeV have been used for the masses of the charm and strange 
quarks respectively. 

The method outlined in Section III allows us to calculate the distortion corrections, which 
comprise the most important contribution to the i^-factor at energies where the quark masses 
are not negligible. At this time we do not interpolate the i^-factor to high-energy results 
involving radiative corrections; the finite temperature radiative corrections derived for the 
quark-gluon plasma |43"1 , [44| may be useful for this purpose. 

For the quark-gluon plasma with a flavor number Nf and N c = 3, lowest-order pertur- 



bative QCD gives a Debye screening length of [38 



^ PQC ^ WtW (24) 
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For a coupling constant a s = g 2 /^ = 0.3 and Nf = 3, the Debye screening length at a 
temperature of 200 MeV is X D ~ 0.4 fm. We examine the cases of X D = 0.2 and 0.4 fm, 
corresponding respectively to temperatures of 400 MeV and 200 MeV in this perturbative 
QCD estimate. We will calculate the associated -fT-factors for charm and strange pairs. 

We show in Fig. 7a the corrective .fT-factor for cc in a color-singlet state for which 
Cf = 4/3 and K(£,r]) > 1. For the c-c interaction, we use the running coupling constant 



2 \ 127T 

as[Q ] = (33-2iV / )ln(gVA 2 )' (25) 
and we assume four flavors and Aqcd = 300 MeV. For X D = 0.4 fm, the iT-factor is much 
greater than unity near the threshold of 3 GeV, and it decreases rapidly as a function of 
energy. The effective coupling constant is approximately a e g = 0.4 in this energy region, 
which gives a Bohr radius of a B »s 0.5 fm. Thus, the screening length of X D = 0.4 fm 
corresponds to i] « 0.8, which is near rj = 0.835 of peak i^-factors, as we can infer from Fig. 
3. The K-f actor oscillates about the Gamow factor (the no-screening limit) which is shown 
as the dashed curve. It is much greater than unity near the threshold. For X D = 0.2 fm the 
if-factor is not as large. This corresponds to a value of rj w 0.4 for this energy region, far 
from rj = 0.835 where the peak of the i^-factor is located. 

We show in Fig. 7b the corrective A'-factor for cc in the color-octet state, for which 
Cf — —1/6 and K(^,rj) < 1. It decreases with increasing X D and varies relatively smoothly 
with energy. The Gamow factor is also shown for comparison. The color-octet corrective 
factor is greater than the Gamow factor, and its deviation from unity is less than the 
deviation of the Gamow factor from unity. 

Because the magnitude of the octet repulsive potential is less than the magnitude of the 
singlet attractive potential, the deviation of the i^-factor from unity is greater for color- 
singlets than for color-octets. 

We have extended our study to strange pair production or annihilation, as strange pair 
production is an important process in high-energy heavy-ion collisions which can be used as 
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a signature of a QGP. To study strangeness production in the quark-gluon plasma, we use 
the running coupling constant proposed by Godfrey and Isgur Q , which is appropriate for 
the energy range near the strange quark pair threshold: 

a s (Q 2 ) = 0.25 exp(-Q 2 ) + 0.15 exp(-Q 2 /!0) + 0.2 exp(-Q 2 / 1 000), (26) 

where Q is in GeV. We display in Figures 8(a) and 8(6) the ss corrective i^-factor in the 
presence of screening for color-singlets and color-octets respectively. 

An interesting difference is observed in the color-singlet K-f actors for charm and strange 
pairs near threshold for the same values for the screening length. For a charm quark pair, the 
corrective i^-factor is larger than for a strange pair, and increases with decreasing energy. 
On the other hand, the color-singlet i^-factor for strange quark pairs decreases slowly with 
decreasing energies. Both systems are studied at similar velocities as we are considering 
regions in which quark masses are important. Although the running coupling constants are 
different because the two systems sample different energy regions at threshold, this is not 
sufficient to account for the difference. 

The difference in K(£, rj) for the charm and strange flavors in a similar velocity ranges 
can be understood in the light of the analysis in Section III. For screening lengths of 0.2 fm 
and 0.4 fm considered here, the values of rj in the cc system are r\ ~ 0.4 and 0.8 respectively. 
The corresponding K(£, rf) have been discussed. For the ss system on the other hand, these 
values of A D correspond to much smaller values of 77, due to the smaller mass of the strange 
quark. We find 77 ~ 0.04 and 0.08 for the two ss cases, and the deviation of K(£,rj) from 
unity is small. This is in accordance with the physical expectation that the interaction is 
suppressed when the range of the potential is small compared to the natural length scale of 
the system a B , so the ss system responds almost as an unperturbed one. 

V. STRANGE AND CHARM QUARK- ANTIQUARK PRODUCTION 

We now consider the effects of screening on cc and ss production cross sections near and 
above the threshold. The basic cross section for qq — > g* — > QQ, (where Q refers to s or c 
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quarks), averaged over initial and summed over final colors and spins can be written as |46 



^(M QQ ) = (l + ^)^l-Vm, (27) 

where r] m = AmQ/MgQ, ijiq is the mass of the quark and Mqq is the invariant mass of the 
produced QQ pair. 

The corresponding expression for the gluon fusion mode, averaged over initial gluon types 
and polarizations and summed over final colors and spins is |46|j2]1 

*-<"<»> = { (1 + "™ + Kl^!) - (I + S"-) } ■ (28) 

To obtain the overall QQ (ss and cc) production, the above cross sections must be convoluted 
with the appropriate quark and gluon distributions. For charm and strange pair production 
in the quark-gluon plasma, these tree-level cross sections must be corrected to take into 
account the effect of Debye screening. 

For QQ production by quark-antiquark annihilation through a virtual gluon mode at 
lowest order, the produced QQ pair is in a color-octet state. The corrective factor to be 
used is the color-octet corrective factor, K qq - = fT(color octet) for the quark-gluon-plasma. 
The corrected cross-section is then 

In QQ production by gluon fusion, QQ pairs are produced with a relative color-octet 
and color-singlet population given by |25|| 

(color - octet) (d* bc /V2) 2 _ 5 

(color - singlet) ~~ (S ah /V3) 2 ~2 ^ ' 

Taking the weights of the population into account, we can write the corrective factor 
for the gluon fusion mode as K gg = [5X(octet) + 2 K (singlet) ]/7. The cross-section for the 
gluon fusion production of QQ pairs, corrected for the effect of screening is then 

W*> = *»3ifc{(l + *. + rfti - (I + V^) } ■ PU 



QQ 
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In Figures 9 and 10 we present the cross sections for cc and ss production corrected for 
screening effects in the quark-gluon plasma, together with the uncorrected cross sections 
for comparison. The cross sections for the case with the color-Coulomb interaction without 
screening are also shown. For ss and cc production by qq annihilation, the intermediate 
virtual gluon selects color-octet states only. The color interaction is therefore repulsive and 
suppresses the cross section. With typical screening lengths of 0.2 and 0.4 fm, the corrected 
cross sections are near the no-screening color-Coulomb values for cc, while for ss production 
they are close to the tree level lowest-order values in accordance with the different values of 
i] for the two flavors. 

The gluon fusion modes on the other hand produce qq pairs in a combination of color- 
singlet and color-octet states, and the larger magnitude of the color-singlet i^-factor dom- 
inates the correction even though the singlet weight factor is lower. Here again, the cross 
sections for cc production approach the no-screening color-Coulomb limit, whereas the cross 
sections for ss production remain close to the lowest-order values, due to the differences 
in T). It may also be remarked that cc production cross sections through the gluon fusion 
mode for \ D = 0.4 fm exhibit a marked rise at small velocities caused by the maximum 
in the K-factor occurring as the lowest bound state emerges into the continuum at r] = 
0.835. On the other hand, the ss production cross section is not appreciably affected by the 
corrective factors in the presence of screened fields. Distortion corrections in such a field are 
not expected to influence ss cross sections in the plasma noticeably. 

VI. CONCLUSIONS AND DISCUSSIONS 

In reactions involving the production or the annihilation of a quark and an antiquark, the 
quark and the antiquark are subject to initial-state or final-state interactions due to their 
mutual interactions. The lowest-order cross sections for these processes can be corrected 
by using an approximate corrective .fT-factor to take into account the q-q interaction. In 
a deconfined quark-gluon plasma, the interaction between the quark and the antiquark is 
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screened. As a consequence, the corresponding ^-factor will be modified. We have studied 
the effect of screening on the corrective if-factor. Our knowledge of such an effect can 
be used to search for the quark-gluon plasma. If the corrective if-factor can be inferred 
from experimental data, then we can use the relation between the color screening and the 
K-factor to deduce the conditions of screening and the possible presence of the deconfined 
phase of strongly-interacting matter. 

In the present investigation we have limited our attention to energies near the q-q thresh- 
old, where the dominant contribution to the corrective if-factor comes from the distortion 
of the wave function by the q-q interaction which is modified by screening. The if-factors for 
the case with screening have been determined by numerical integration of the non-relativistic 
Schrodinger equation for q-q with a color- Yukawa potential. The if-factor gives an enhance- 
ment for color-singlet qq states, due to the attractive nature of the color-singlet potential, 
but a suppression for color-octet qq states due to the repulsive nature of the color-octet 
interaction. The color-singlet enhancement is larger in magnitude than the color-octet sup- 
pression because of the difference in the two strengths specified by the color factors of 4/3 
and —1/6 respectively. 

The corrective if-factor for the Yukawa interaction, — a cS e~ r / XD /r, depends on two di- 
mensionless parameters: the usual 'Coulomb' parameter £ = a e g/v, and the screening length 
parameter rj = \ D /a B . We have evaluated K(£,r}) for a large range of of £ and 77, to study 
how the if-factor depends on the q-q relative velocity and the screening length. The function 
K(£,r)) expressed in the dimensionless parameters £ and rj can be used for other reactants 
interacting through a Yukawa potential. 

For attractive Yukawa potentials we observe prominent peaks of the If -factor as a func- 
tion of the screening length parameter 77. The peaks are located at r] = 0.835 and 7] = 3.23, 
corresponding to lowest two s-wave qq bound states emerging into the continuum to become 
qq resonances as the screening length decreases. These qq resonances can be called "screen- 
ing resonances" as they arise from the screening of the q-q interaction in the quark-gluon 
plasma. We have shown that their behavior is compatible with resonances expected for 
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potentials with finite range. This is also in accord with results obtained for the scattering of 



an electron from an atom p5|j36|j, and is similar to the manifestation of resonances through 



the s-wave strength functions in low-energy neutron scattering from nuclei ||40|| . The explicit 
location of the peak at X D = 0.835a s coincides with the screening length at which the lowest 
qq state becomes unbound ||. 

We have calculated the K-factors for two typical choices of the Debye screening length 
(0.2 fm and 0.4 fm), corresponding to plasma temperatures of approximately 400 and 200 
MeV respectively. They are given separately for color-singlet and color-octet q-q states and 
can be used to correct the lowest-order cross sections. We have used these i^-factors to 
investigate the effects of screening on the cross sections for charm and strange pair produc- 
tion. While the corrections to the cross section in the no-screening color-Coulomb limit are 
of similar magnitude for ss and cc pairs in the same velocity range, they are considerably 
different for the two systems in the presence of screening with a color- Yukawa interaction. 
This arises because the Bohr radius is much smaller for the c-c system than for the s-s 
system. A screening length of 0.4 fm would correspond to rj rs 0.8 for c-c but only 0.08 for 
s-s. The screening length paramter 77 = 0.8 happens also to be near the location where a c-c 
resonance occurs at E m 0. Thus, c-c color singlet production is much enhanced for A = 0.4 
fm. There is no such enhancement for s-s production at this screening length. 

Since the single virtual gluon in qq annihilation leads to color octet pair production, 
the lowest order rates for qq — > cc(ss) are larger than the corrected ones because the octet 
correction is suppressive. On the other hand, qq production by gluon fusion can proceed 
through color-singlet or octet states, the weight being given by 2/5 for singlet to octet. 
Despite the lower weight, the larger magnitude of the color-singlet i^-factor results in the 
attractive potential dominating, so that the net effect of the correction is an enhancement 
of gluon fusion production. 

It is also noteworthy that the peak of the fT-factor in the cc system corresponds to the 
unbinding of the lowest bound state and its emergence as a cc resonance. It is manifested 
as a rise in the gluon fusion cross sections near threshold, for a screening length of 0.4 fm, 
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which is near the screening length required for a peak in the i^-factor, as shown in Fig. 
9. Gluon fusion being the dominant mode for heavy-quark production, the combined effect 
is to have a peak in the cc production cross section just above the cc threshold, with a 
relatively narrow width in energy. Thus, the occurrence of a cc resonance will be signalled 
by a large enhancement of cc production just above the threshold. This is in contrast to ss 
production for X D = 0.2 — 0.4 fm, for which the corresponding screening length parameters 
of rj = 0.04 — 0.08 are very small and are far from r\ = 0.835 for a qq resonance. 

From our results, the occurrence of qq screening resonances may provide a way to search 
for the quark-gluon plasma. The qq resonances give rise to prominent peaks of the A'-factor 
as a function of the screening length parameter, which is the ratio of the screening length 
X D to the Bohr radius a B . The screening length is inversely proportional to the temperature 
||38|| . Thus, qq resonances lead to prominent peaks of the i^-factor as a function of the 
plasma temperature. We have seen in Fig. 9 that large values of the X-factor near the 
threshold give rise a narrow peak in the heavy-quark production cross section just above the 
threshold. The occurrence of a qq resonance will be accompanied by a much enhanced qq 
production cross section just above the threshold, and the enhancement will be a function 
of the temperature. 

The search for qq screening resonances in the quark-gluon plasma can make use of the 
peaks of the i^-factor at rj = 0.835 and rj =3.23. The resonance at rj = 3.23 may not lead to 
realizable enhancements because it corresponds to temperatures much below the estimated 
quark-gluon plasma transition temperature (of approximately 150 — 200 MeV). Using the 
perturbative QCD estimates, the screening length parameter of rj = 0.835 corresponds to a 
c-c resonance at T cE pa 165 MeV and a bb resonance at T b i pa 393 MeV. These T c5 and T b i 
estimates from PQCD are approximate and uncertain, as lattice gauge theory gives Debye 
screening lengths of about half of PQCD estimates. (See Fig. 7 of ^7|.) The Debye screen- 
ing length needs to be determined by experimental searches for these cc and bb resonances 
using the peaks in the iT-factors. The PQCD estimates are useful only as a rough guide. 
Thus, for T pa T c5 (approximately 165 MeV), and if the quark-gluon plasma is produced, 
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r)(cc) = 0.835 and there will be a large enhancement of cc production near the threshold, 
much greater than what one expects from the lowest-order cross sections (Fig. 3). At tem- 
peratures far from the cc resonance temperature T cS there is no such large enhancement of 
cc production near threshold. At T m T b i (approximately 393 MeV), if the quark-gluon 
plasma is produced, rj(bb) = 0.835 and there will be a large enhancement of bb production 
near the threshold, again much greater than expected from the lowest-order cross sections. 
At temperatures far from the bb resonance temperature T b i, there is no such large enhance- 
ment of bb production near the thresholds. Temperature dependence of this type arises 
from the nature of screening between the interacting heavy quark and its antiquark partner, 
which is an important property to identify the deconfmed quark-gluon plasma. A search 
for the quark-gluon plasma using heavy-quark resonances will require the measurement of 
the production yield of heavy-quark pairs near the threshold, and a method to estimate the 
temperature of the environment in which the production takes place. The enhancement will 
occur either for the production of heavy-quark pairs by the collision of the constituents of 
the thermalized quark-gluon plasma, or by the collision of the partons in nucleon-nucleon 
collisions in a quark-gluon plasma environment. 

It is interesting to note the recent report from the NA38 and the Helios-3 Collaborations 
of an excess of dileptons at a dilepton invariant mass of 1.5 to 2.5 GeV [^8-50|. The NA38 



Collaboration suggests that this may be due to an excess production of charm pairs in 



high-energy nucleus-nucleus collisions |[49|| . The estimate from our work gives cc screening 
resonance occurring at a temperature T c5 about 165 MeV, which is within the experimental 
range of temperatures as encountered in the NA38 and the Helios-3 experiments. It will be 
of great interest to have a direct measurement of charm pair production near the threshold 
for the reactions studied by NA38 and Helios-3, to see whether there are indeed excess charm 
pairs near the threshold, and whether these excess charm pairs arise as a manifestation of 
a cc resonance emerging just above the threshold in a quark-gluon plasma, or from some 
other processes (TH . 



After the work in the present paper was completed, a recent study of the effect of a re- 
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pulsive screened Coulomb interaction between two pions in pion interferometry was brought 



to our attention [^J. The modification of the Gamow factor due to screening was obtained 
in |HJ by using a quasi-classical approximation to barrier penetration, whereas the results 
obtained here (Section III) for repulsive screening potential follows from exact integration of 
the Schrodinger equation. For this special case of a repulsive screened potential, our exact 



results complement the quasi-classical approximate solutions of 51 
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FIGURES 

FIG. 1. The corrective factor if (£,77) plotted as a function of £ for different values of the 
dimensionless screening parameter 77. The parameter £ is positive for an attractive Yukawa potential 
and is negative for a repulsive Yukawa potential. The Gamow factor for the case of Coulomb (or 
color-Coulomb) interaction is also shown as the dashed curve. 

FIG. 2. The corrective factor if (£, rj) plotted as a function of £ for 1 < £ < 3 and various values 
of 77. The Gamow factor for the case of Coulomb interaction is also shown as the dashed curve. 

FIG. 3. (a) The corrective factor if(£, rj) plotted as a function of 77 for 0.4 < 77 < 2 and different 
values of £, (b) the ratio -fT(£,7/)/(Gamow factor) as a function 77 for different values of £. 

FIG. 4. The phase shift 5q as a function of energy for 77 near 0.835. 

FIG. 5. (a) The corrective factor K (£,77) plotted as a function of 77 for 1.5 < 77 < 5.0 and 
different values of £, (b) the ratio if(£,77)/(Gamow factor) as a function 77 for different values of £. 

FIG. 6. (a) The corrective factor if (£, rj) plotted as a function of £ for — 3 < £ < and different 
values of 77. Negative values of £ represent a repulsive Yukawa interaction. The Gamow factor 
for the case of Coulomb interaction is shown as the dashed curve, (b) the ratio if(£, rj) / (Gamow 
factor) as a function £ for different values of 77. 

FIG. 7. The corrective factor if(£, 77) for a cc pair with screening length A D =0.2 and 0.4 fm, 
as a function of the center-of-mass energy y/s. The Gamow factor is shown as a dashed curve for 
comparison. Fig. 7(a) is for color-singlet cc states, and Fig. 7(b) is for color-octet cc states. 

FIG. 8. The corrective factor if (£,77) for an ss pair with screening length A D =0.2 and 0.4 fm, 
as a function of the center-of-mass energy y/s. The Gamow factor is shown as a dashed curve for 
comparison. Fig. 8(a) is for color-singlet ss states, and Fig. 8(b) is for color-octet ss states. 
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FIG. 9. Cross sections for cc production (a) by gluon fusion, and (b) by qq annihilation, with 
screening lengths of X D = 0.2 and 0.4 fm. The lowest-order cross section and the cross section in 
the case of color-Coulomb c-c interaction are also shown. 

FIG. 10. Cross sections for ss production (a) by gluon fusion, and (b) by qq annihilation, with 
screening lengths of \ D = 0.2 and 0.4 fm. The lowest-order cross section and the cross section in 
the case of color-Coulomb s-s interactions are also shown. 
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